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Abstract. Experimental and numerical studies of the velocity field of dark solitons 
in Bose-Einstein condensates are presented. The formation process after phase 
imprinting as well as the propagation of the emerging soliton are investigated using 
spatially resolved Bragg-spectroscopy of soliton states in Bose-Einstein condensates of 
^^Rubidium. A comparison of experimental data to results from numerical simulations 
of the Gross-Pitaevskii equation clearly identifies the flux underlying a dark soliton 
propagating in a Bose-Einstein condensate. The results allow further optimization 
of the phase imprinting method for creating collective exitations of Bose-Einstein 
condensates. 
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1. Introduction 

With tlie realization of Bose-Einstein condensation (BEC) in dilute atomic gases PIEIISI, 
a well controllable and easily observable macroscopic wavefunction has become accessible 
for a wide variety of fundamental experiments. Since the early experiments on 
oscillations of the condensate |31 Ej and phonons [Hj, excitations of the condensate 
wavefunction have been subject to intense investigation. Signatures of superfluidity 
in Bose-Einstein condensates were first observed in the suppression of excitations below 
a critical excitation velocity [3 |H] and in the analysis of scissors mode excitations [H]. 
Recently another striking aspect of superfluidity in BEC has been studied by exciting 
vortices in two-component condensates ^lOj and vortices and vortex arrays in stirred 
condensates [TTl [T21 [13 Cll 113 CHI • Vortices as topological excitations of the condensate 
wavefunction represent particularly interesting phenomena and are still a very active 
area of research. Dark solitons in Bose-Einstein condensates are in several respects the 
one-dimensional counterpart to vortices. While the phase of the vortex wavefunction 
winds up around a zero density 'core', in a dark soliton state regions of constant phase 
are separated by a minimum in density and a steep phase gradient. Dark solitons 
are fundamental excitations in one dimensional systems described by nonlinear wave 
equations and thus a solution of the Gross Pitaevskii equation, which describes the 
condensate wave function In Bose-Einstein condensates, the creation of dark 

solitons using a phase imprinting method jtlSj has been investigated in a strongly 
elongated geometry in BECs of ^^Rb ^Hl and in nearly spherical BECs of ^^Na |2()j . 
Recently bright solitons in Bose-Einstein condensates with attractive interactions 
|2H 1221 ^'^'^ filled solitons j2Sl have also been observed . 

The experiments on dark solitons have so far concentrated on the temporal evolution 
of the density distribution of the Bose-Einstein condensate, which showed the motion 
and decay of these excitations. In this article we extend these investigations to the 
velocity field which underlies a dark soliton state in a Bose-Einstein condensate. 

2. Properties of dark solitons 

In general Solitons are fundamental wave-packet-like excitations of one dimensional 
nonlinear systems. One characteristic of them is that they do not disperse but preserve 
their shape due to a balance between dispersive and stabilizing effects in these systems. 
Dark solitons in Bose-Einstein condensates form density minima, which self stabilize 
due to a balance of the repulsive mean field interaction which would be minimized by 
compressing the minimum and quantum pressure. The quantum pressure is related to 
the kinetic energy associated with the dark soliton velocity field, which is minimized by 
widening the minimum. 

More formally a dark soliton in a homogeneous Bose-Einstein condensate of density 
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Figure 1. Schematic of the density distribution (a), phase distribution (b), and 
the resulting velocity distribution (c) of a dark soliton propagating in a BEC. 
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with the position and velocity of the dark soliton, the correlation length /q = 
{Anario)'^/'^ , and the speed of sound Cg = y/Ananoh/m, where m is the atom mass. 
Note that a so called black soliton with zero density at the minimum and a vr phase 
jump is obtained for Vk = 0, while for all other cases the condensate wave function has 
nonzero amplitude at the position of the density minimum which now moves along the 
condensate and is associated with a finite phase gradient. 

Consequently dark solitons in Bose-Einstein condensates are density minima 
associated with a steep phase gradient of the condensate wavefunction. This gradient of 
the phase distribution <l>(r ) of the macroscopic Bose-Einstein condensate wavefunction 
defines a probability current density, which can be interpreted as a superfluid velocity 
field f (r ) = ^V$(r ). Thus, the steep gradient of the phase at the dark soliton leads 
to a localized peak in the velocity distribution of the BEC. The corresponding BEC 
density, phase and velocity distributions are shown schematically in Fig. 

Equation ^ is strictly valid only for infinite homogeneous one- dimensional systems 
but many qualitative features remain valid in more realistic experimental cases. However 
the effect of the trapping potential, the formation process and the dynamics during time- 
of-flight expansion are not included. In order to model these cases more accurately, 
we numerically solve the Gross-Pitaevski equation for the case of inhomogenous 
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condensates, trapped in static potentials. Deviations from the ID behaviour in elongated 
and nearly spherical traps are treated e.g. in [23 123 12Z1 ■ 

3. Creation of dark solitons 

For the experimental investigation of the dark soliton velocity field we first produce 
ground state Bose- Einstein condensates of ®^Rb as described in [28j. The condensates 
typically contain 1.5 x 10^ atoms in the {F=2, mp=+2) state, with less than 10% of the 
atoms being in the thermal cloud. The fundamental frequencies of our static magnetic 
trap are = 2it x 14 Hz and uj± = 271 x 425 Hz along the axial and radial directions. 
The condensates are cigar-shaped with the long axis (x-axis) oriented horizontally and 
an aspect ratio of A i=5i 30, providing an essentially one dimensional potential for the 
evolution of the dark soliton on the time scale of the experiment. 

In our experimental setup we create dark solitons using a phase imprinting 
method [IHl dl. The phase distribution of the condensate is changed by using the 
dipole potential of a far detuned light field (A =532 nm, I ^ 20W/mm^), illuminating 
one half of the condensate for a time tp on the order of 20 /is. Note that this time has 
to be shorter than the correlation time of the condensate, tc = h/fi, in order to obtain 
pure phase imprinting, where /i is the chemical potential and in our case t^ ~ 55 /is. 
This procedure creates a moving dark soliton and a counterpropagating density wave, 
which become well separated after some time, tev, of further nonlinear evolution in the 
trapped Bose-Einstein condensate, as shown in Fig. |21a). 

Figure |21 b) shows a typical density-distribution of a dark sohton state created by 
phase imprinting after some evolution time in a trapped BEC obtained by the numerical 
solution of the GPE in ID using the split operator method This simulation contains 
a finite optical resolution for the phase imprinting field, leading to a potential gradient 
with a width of 2 fim at the edge of the illuminated realized experimentally. 

Due to this gradient a density wave and some higher order excitations moving to the 
right are created in addition to the dark soliton. The dark soliton and the density wave 
move in opposite directions, however the dark soliton is much narrower due to its self 
stabilization in the nonlinear medium. Note that the velocity field in the region of the 
dark soliton has only components of negative velocities, v{x) < for a positive soliton 
velocity Vk > 0. This can be intuitively understood by the fact that a movement of the 
density minimum necessitates a transfer of matter in the opposite direction. 

4. Bragg spectroscopy 

Bragg diffraction has been used as a tool to manipulate Bose-Einstein condensates jSH] 
in a number of experiments. It was used to evaluate the velocity distribution of a 
BEC within the magnetic trap and during time-of- flight experiments jSUj. It has also 
served as a beam-splitter in an interferometer [HI], for interferometric studies of the 
phase evolution of BECs jH21 and even as an outcoupling mechanism for an atom laser 
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Figure 2. a) Absorption images of Bose-Einstein condensates which evolved 
in the magnetic trap for a time, tg^ after phase imprinting. The images were 
taken after an additional time-of-flight of 4 ms, leading to density minima for 
the dark soliton as well as for a density wave JH]- b) Numerical simulation of 
the density distribution of a dark soliton propagating in a BEC. The density 
data shown was obtained for the same parameters as used for the velocity field 
in Fig. ini 

PH] . In addition Bragg diffraction was employed as a tool for the realization of exciting 
experiments on non-linear effects in BECs jSlj and their possible application to create 
entanglement jHS]. However, a complete review of the use of Bragg spectroscopy in 
experiments on BEC is beyond the scope of this paper. 

Fig. 121 schematically shows the Bragg diffraction process used in the experiments 
described here (a more detailed explanation can be found in j23|)- Two 
counterpropagating laser beams induce a two photon transition between two momentum 
states which differ by two photon momenta, 2hk. Thus the atoms which have been in 
resonance with the Bragg-transition will be separated from the remaining atoms by 
Ax = , where ttof is the time-of-flight after the Bragg pulse. Since the resonance 
conditon depends on the initial velocity of the atom, a given difference frequency of 
the laser beams addresses a specific velocity class within the condensate. The width of 
this resonance corresponds to the Fourier transform of the applied laser pulse envelope. 
For ^''Rb and counterpropagating laser beams the velocity dependence of the resonance 
condition is given by ~ 0.39^^^. 

Experiments employing Bragg diffraction on BEC can be divided in two general 
areas. By using short, intense pulses all velocity classes within a BEC can be accessed 
and transferred into a moving state. This process corresponds to a variable, coherent 
beamsplitter for matterwaves. On the other hand, long pulses are selective to the 
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Figure 3. a) Sciiematic of position resolved Bragg-spectroscopy of a dark soliton 
state in a BEC. b) Energy level scheme for the Bragg transition between ground 
state levels with velocities v and v + where k is the wave number of the 

m 

laser field. The resonance condition for a given frequency difference 5 is only 
fulfilled for a certain velocity v. 

velocity in the BEC and distinct velocity classes can be addressed. This process is used 
to examine the flux underlying dark solitons in Bose-Einstein condensates. We show 
that it is also possible to extract spatially resolved information from these measurements, 
rather than only counting the outcoupled number of atoms. 

The experimental sequence starts with the creation of a Bose-Einstein condensate 
in the magnetic trap followed by phase imprinting for the creation of the soliton (as 
discussed above). After phase imprinting the condensate evolves in the magnetic trap 
for a time tev which is varied in order to observe the soliton excitation during formation 
and later propagation. The Bragg scattering takes place after switching off the magnetic 
trapping potential. For the experiments a free expansion time, t^xp = 2 ms, is used after 
release from the trap in order to reduce the atomic density before applying the Bragg 
pulse. Collisions between atoms in different Bragg orders are thus negligible. In order 
to suppress spontaneous processes the frequencies of the Bragg lasers are detuned from 
the transition |F = 2) \F') by 6.8 GHz. The Bragg pulse is usually applied for a 
duration of tBragg = 250 yus, and a time-of-flight of ^tof = 13 ms after the pulse is used 
before detection (tBragg = 500 /is and troF = 11 ms for the BEC spectroscopy in Fig. 
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El). 

5. Velocity selective spectroscopy 

The velocity selective Bragg scattering process described above is a widely applicable 
tool to investigate the behaviour of Bose- Einstein condensates. In the following we first 
observe the formation process of dark solitons due to the phase imprinting method. This 
is achieved by monitoring the spatial velocity almost immediately after the phase has 
been imprinted. In a second set of experiments the soliton is observed after a longer 
evolution time, tev, in the magnetic trap, in order to investigate the associated matter 
flux. 

5.1. Soliton formation 

Fig.m shows absorption images of atomic clouds after the above Bragg spectroscopy 
scheme has been applied with tev = 0.2 ms and different resonance conditions. The 
upper image (i) was obtained for an off- resonant Bragg pulse and thus just corresponds 
to a time-of-fiight image of the BEC after phase imprinting. Image (ii) shows the case 
of a Bragg pulse resonant with the resting atoms, which thus obtain a velocity, 2hk, and 
appear to the right of their position in image (i) due to their horizontal movement during 
the time-of-fiight period of 13 ms after phase imprinting. The moving atoms within the 
cloud are not resonant with the Bragg pulse and remain on the left. They correspond 
to the density wave and dark soliton. The velocity of these atoms can be determined 
to be negative, as expected, by scanning the Bragg frequency difference. The density 
distributions corresponding to Bragg scattering resonant with negative velocity atoms 
are shown in image (iii), with the resting atoms on the left and the moving atoms on 
the right. 

The simulations show that the velocity field in the Bose-Einstein condensate after 
the short evolution time, tev = 0.2 ms corresponds almost exclusively to the velocity 
field induced by the phase imprinting pulse. Fig.El shows calculated density and phase 
distributions in the magnetic trap for various times after the phase imprinting pulse. 
At tev = there is no density change but only a phase gradient in the center of the 
condensate, which corresponds to a velocity field. The associated flow quickly leads to 
the formation of a density maximum, moving to the left and a density minimum moving 
to the right. These two start to separate after approximately tev = 0.5 ms with a 
subsequent dispersive spreading of the density maximum and a narrowing of the density 
minimum which leads to the dark soliton state. Note that the initial phase gradient is 
split between the density wave and the dark soliton since both are associated with some 
flow of matter. 

The evolution of the wavefunction in Fig.El clearly shows that the mechanism of 
phase imprinting with finite resolution relies on the joint action of a phase shift and a 
local removal of atoms due to a momentum transfer. The momentum imparted on some 
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Figure 4. a) Absorption images for different detunings 5 of the Bragg-pulses: 
(i) off-resonant, (ii) resonant with velocities v = 0, (iii) resonant with velocities 
V < 0. These images were obtained for a short evolution time in the magnetic 
trap of tev = 0.2 ms. 
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Figure 5. Density and phase of the condensate wavefunction from numerical 
solutions of the Gross Pitaevskii equation for short evolution times in the 
magnetic trap after phase imprinting. For these simulations the width of 
the light field edge was assumed to be 2 ^um. These simulations show that 
for the experimental parameters the soliton and density wave are not yet 
distinguishable. Note that in this case the simulation only shows the central 
part of the condensate and doesn't include any time of flight. 




Figure 6. Left: velocity distribution of a dark soliton state after phase 
imprinting with a phase difference At/; = tt, an evolution time in the magnetic 
trap tey = 3.5 ms, and tj-oF = 2 ms of time of flight from a numerical simulation 
of the Gross-Pitaevskii equation. Right: corresponding Bragg spectroscopy 
images for different resonance conditions v = +0.5mm/s, v = — 0.7mm/s 
and V = — l.lmm/s (from top to bottom). The height of the vertical bars 
indicates the width of the Bragg resonance, thus representing the velocity 
classes, which will be scattered during the Bragg pulse. 



parts of the wave function, assists in the later formation of a dark soliton due to the 
removal of density from the region of the soliton. The images in Fig.E] confirm these 
numerical results. The imparted momentum can be explained classically by the finite 
width of the edge of the phase imprinting field, which produces a classical force due to 
the gradient of the light induced AC-Stark shift potential. 

5.2. Flux associated with a dark soliton 

At short evolution times the dark soliton is masked by other excitations induced by 
the phase imprinting pulse, but after a few ms it becomes apparent due to its self 
stabilization and slow propagation. An evolution time of t^v = 3.5 ms is sufficient 
to fully seperate the soliton and the density wave as shown by the Bragg scattering 
experiments and the numerical simulation in Fig. IHl 

From the numerical simulation one expects to find local negative velocity 
components around v ~ — 0.5m/s at the positions of the density wave as well as of 
the dark soliton. The top right image shows nearly resting atoms scattered to the 
right, while two distinct lines of moving atoms remain on the left: the atoms associated 
with the dark soliton and the density wave. This scattering behaviour is inverted for a 
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Figure 7. Vertically summed cuts through Bragg spectroscopy images of a 
Bose-Einstein condensate without phase imprinting. The Bragg pulse duration 
for this data was 500 ^s. The Bragg resonance frequency is changed by 200 Hz 
(corresponding to Av ^ 0.08 mm/s) between subsequent cuts. (Note that for 
technical reasons the position of the scattered atoms is on the other side as for 
the other images.) 

negative velocity Bragg resonance, now scattering the two lines to the right and leaving 
the resting atoms on the left. No atom is effected by a Bragg pulse resonant with larger 
negative velocities in the bottom image, showing the limited range of the velocity field 
in the matter wave. This is clear evidence of the matter wave flux associated with a 
moving dark soliton. 

6. Mapping the velocity field 

In addition to the above examples of Bragg scattering for a few distinct resonance 
velocities, which clearly show the signature of dark soliton excitations, we use this 
method to fully map the underlying velocity field. 

We first explain and verify the scheme by mapping the spatial velocity distribution 
of a pure Bose-Einstein condensate before applying it to a dark soliton excitation. 
The starting point is a systematic series of Bragg scattering data of a Bose-Einstein 
condensate after 2 ms of time of flight, as shown in Fig. [7| By taking a vertical cut 
through the data sets in the region of the scattered peak one obtains the amount of 
scattered atoms at a fixed position as a function of Bragg detuning. The maximum 
of this distribution then directly leads to the mean velocity of the atoms scattered to 
this position. Their original position in the condensate is then retrieved by accounting 
for their time-of-fiight displacement due to the imparted momentum and their original 
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Figure 8. Velocity field of a freely falling Bose-Einstein condensate as 
determined from the series of Bragg scattering data in Fig. 

velocity. Repeating this procedure for the position range of the scattered peak this 
procedure directly gives the velocity field of the condensate wavefunction. In addition 
to analyzing the scattered peak, we cross check the velocity field data by generating the 
same data from the inverse behaviour of the non-scattered cloud. 

Fig. ISlshows the extracted velocity as a function of position in the Bose condensate. 
The result is fitted with a linear velocity distribution according to the expected 
behaviour in the Thomas- Fermi limit. The good agreement confirms the validity of 
this analysis. 

By applying the above method to the dark sohton excitation discussed in section 
15.21 we find the underlying velocity field as shown in Fig. IHl 

The comparison of the experimental data to the theoretical curve shows 
good agreement within the finite experimental resolution, which demonstrates the 
applicability of position resolved Bragg spectroscopy to the investigation of excitations 
in Bose-Einstein condensates. 

7. Conclusions 

In conclusion we have investigated the matter wave dynamics due to phase imprinting 
and the matter wave flux underlying a moving dark soliton. Both experiments and 
theoretical predictions clearly show that the dark soliton movement is associated with 
a particle flux in the opposite direction. It was demonstrated that Bragg diffraction 
methods can be used to analyse the spatial phase and the density distribution of dark 
soliton states in Bose-Einstein condensates. 
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Figure 9. Velocity field underlying the dark soliton excitation in section 15.21 
in comparison to the numerical simulation from Fig. El The fit through the 
experimental data is a spline fit to guide the eye. The negative velocity peaks 
(a) and (b) correspond to the density wave and the dark soliton respectively. 
Note that the averaging effect of the experimental data is mainly due to the 
finite optical resolution of « 15 //m used for this set of measurements. 

Future improvements of the tool of the position resolved Bragg spectroscopy 
technique rely on the use of longer pulses and higher resolution in the detection system. 
This should allow for precision measurements on many kinds of excitations in Bose- 
Einstein condensates. In this respect the investigation of the interplay between the 
imprinted phase and the associated imparted momentum promises important advances 
in the creation and study of excitations in Bose-Einstein condensates. 
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